1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.02
## LDevsum 1.01 1.03
## dh0 1.50 2.15
## dh1 1.66 2.63
## dh2 1.41 2.22
## dl0 1.92 3.06
## dl1 1.60 2.49
## dl2 1.39 2.11
##
## Multivariate psrf
##
## 1.83
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1158.7199 | 22017.5632 |
| DIC3 | 1122.8126 | 22085.9396 |
| PWAIC | 26.4105 | 187.0915 |
| WAIC | 1136.2745 | 22098.6314 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.01
## LDevsum 1.00 1.01
## dh0 2.55 8.40
## dh1 9.41 17.85
## dh2 5.53 13.34
## dl0 1.12 1.33
## dl1 1.17 1.47
## dl2 1.05 1.14
##
## Multivariate psrf
##
## 8.43
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1380.41175 | 30292.6188 |
| DIC3 | 1244.84499 | 30382.9877 |
| PWAIC | 90.14513 | 139.2839 |
| WAIC | 1313.83189 | 30387.0876 |